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The overlap hypercube fermion is constructed by inserting a lattice fermion 
with hypercubic couplings into the overlap formula. One obtains an exact 
Ginsparg- Wilson fermion, which is more complicated than the standard over- 
lap fermion, but which has improved practical properties and is of current 
interest for use in numerical simulations. Here we deal with conceptual as- 
pects of the overlap hypercube Dirac operator. Specifically, we evaluate the 
axial anomaly and the index, demonstrating that the correct classical con- 
tinuum limit is recovered. Our derivation is non-perturbative and therefore 
valid in all topological sectors. At the non-perturbative level this result had 
previously only been shown for the standard overlap Dirac operator with 
Wilson kernel. The new techniques which we develop to accomplish this are 
of a general nature and have the potential to be extended to overlap Dirac 
operators with even more general kernels. 
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1 Introduction 



In gauge theories with fermions, the index of a Dirac operator plays an impor- 
tant role. It is given as the difference of the number of left- and right-handed 
zero modes, and due to the Atiyah-Singer Index Theorem ^ it can be iden- 
tified with the topological charge of the gauge field. Therefore the same 
quantity also provides the integrated axial anomaly. 

On the lattice it is a non-trivial question if these quantities can be re- 
covered. In particular it is not clear a priori if we obtain the correct axial 
anomaly in the classical continuum limit. The traditional lattice formulation 
by means of the Wilson Dirac operator D\y does not allow for the Index 
Theorem to be adapted: lattice gauge configurations do not have natural 
topological sectors, and the Wilson fermion does not have exact zero modes. 
Nevertheless the axial anomaly can be reproduced j2] ; a necessary condition 
for this property is the absence of species doublers. 

The situation is different for overlap fermions They have good chiral 
properties according to their origin from the overlap formalism |3115|, which 
are reflected by the fact that their lattice Dirac operator Dqv obeys the 
Ginsparg- Wilson relation (GWR) IHIIZI- In an even dimension 2n it reads ^ 

Dovl2n+l + l2n+lDov = —Dovl2n+lDov ■ (1-1) 

m 

Here m is a parameter which controls topological properties and the number 
of fermion species described by the corresponding lattice fermion action. The 
GWR turns into the standard condition for chiral symmetry in the continuum 
limit. Moreover, even on the lattice an exact chiral symmetry exists jH]- It 
is lattice modified by a local term of 0(a). Explicitly, the variation of the 
spinor fields under the lattice modified chiral transformation can be written 
in the form "H] 

5^) = ^Dov)l2n+l , Sip = 'y2n+l(l - ^Dov)lp , (1-2) 

which leaves the lattice action invariant. 

The Ginsparg- Wilson relation excludes additive mass renormalization. 
The index is well defined, since the exact zero modes have a definite chirality 

^We refer to a Euclidean lattice of spacing a, and the fermions belong to some unitary 
representation of the (unspecified) gauge group. 
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[7j. ^ The question arises if the correct axial anomaly is reproduced in the 
continuum limit. In fact, the existence of some anomaly is obvious since 
the fermionic measure is not invariant under the transformation ()1.2j) . The 
resulting anomaly takes the form 

Aix) = 2iq{x) , (1.3) 

where q{x) is the topological charge density, 

q{x) = -^tY(j2n+lDov){x,x) . (1.4) 

Generally we define the density 0{x,y) of a lattice operator by 

0^|J{x) = a'''J2&ix,y)^P{y) . (1.5) 
y 

Then the operator trace can be expressed as 

Tr0 = a2"^tr0(x,x) , (1.6) 

X 

where tr sums over spinor and gauge indices. In particular the index is given 
by ^ 

mdex{Dov) = -—Tt {-f2n+iDov) = a^"^g(x) . (1.7) 

The question if one obtains the correct expression for the axial anomaly 
has been studied in a number of papers. This issue was first addressed in Ref. 
[0], where a perturbative evaluation of the anomaly was given for a general 
Dirac operator satisfying the Ginsparg- Wilson relation. ^ Furthermore, a 
perturbative evaluation of the axial anomaly for very general lattice Dirac 
operators has been carried out in Refs. jTU]- In all cases the anomaly was 
found to reduce to the correct continuum expression if the operator is local 
and free of species doubling. However, the use of perturbation theory restricts 
the validity of these arguments to the sector of topologically trivial gauge 
fields. In fact, the question for which Ginsparg- Wilson fermions the correct 

^In the case of the overlap Dirac operator the index coincides with the overlap topo- 
logical charge For the fixed point fermion, which also solves the GWR, it coincides 
with the classically perfect topological charge .7^. 

■^This perturbative evaluation was reconsidered in the context of modern developments 
of the Ginsparg- Wilson relation 
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anomaly holds also in topological sectors different from zero seems to be 
non-trivial in the light of Ref. ^1], which presents an example where this 
is apparently not the case. So far, the only Ginsparg- Wilson operator for 
which the correct continuum limits of the axial anomaly and index have 
been established non-perturbatively is the standard overlap Dirac operator, 
which is given by substituting the Wilson-Dirac operator D = D\y into the 
overlap formula 



The continuum limit of the axial anomaly in this case has been studied 
explicitly in Refs. [IS UHl Ell EH] • In particular, a rigorous non-perturbative 
demonstration that the anomaly and the index have the correct continuum 
limit in all topological sectors when m is in the physical (doubler-free) region 
was given in [121 EI . 

Currently there is interest in non-standard overlap Dirac operators ob- 
tained by inserting more general lattice Dirac operators D into the overlap 
formula (|1.8|) . The background and motivation for this is discussed further 
below. In particular, overlap Dirac operators where the input D is a hy- 
percubic fermion (HF) operator have been the focus of attention.^ In the 
light of this development it is pertinent to show at the non-perturbative level 
that the axial anomaly and index for these non-standard overlap Dirac op- 
erators also have the correct continuum limit in all topological sectors. As 
mentioned above, this has so far only been shown for the standard overlap 
operator where the input D is the Wilson-Dirac operator. The purpose of 
the present paper is to establish this result at the non-perturbative level 
for non-standard overlap Dirac operators of specific current interest, namely 
those for which the input D in the overlap formula is a HF operator. To do 
this we follow the rigorous non-perturbative approach of Refs. ||13| il4j for 
the standard overlap Dirac case, and take inspiration from Ref. where a 
topological description of the anomaly coefficient as the degree of a certain 
map was derived which greatly facilitates its evaluation. However, the key 
technical parts of the arguments in those papers are specific to the standard 
overlap case — they rely on the explicit form of the Wilson kernel and do not 
have a straightforward generalization to more general kernels. Therefore we 

^HF operators are generalizations of the Wilson-Dirac operator which couple all sites 
within a lattice hypercube. (Recall that the Wilson-Dirac operator couples only nearest 
neighbor sites.) 
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have had to develop new techniques and formulae to handle the more general 
HF case. In fact, our techniques are of a general nature and have the po- 
tential to be used for even more general kernels. (The case of overlap Dirac 
operator with completely general kernel has further complications though, 
and is postponed to a later analysis.) 

Let us now discuss the background and motivation for considering the 
non-standard overlap operators mentioned above. Their use was suggested 
in Ref. ^Hl- The motivation is to improve other properties of the overlap 
Dirac operator — beyond chirality — which are also of importance for a 
lattice fermion formulation, such as the quality of scaling, locality and of 
approximate rotation invariance. We emphasize that chiral properties of the 
overlap operator Dov continue to hold for any input lattice Dirac operator D 
(free of species doubling) in ()1.8|) . since Dov satisfies the GWR for any such 
choice. The basic idea is to construct a short range, doubler-free lattice Dirac 
operator D for the input into the overlap formula which has the following 
properties: (i) good scaling and approximate rotational invariance, and (ii) 
good chirality in the sense that D is an approximate solution of the GWR. 
The property (ii) suggests that the overlap Dirac operator Dov obtained from 
inserting D into the overlap formula will inherit to a large degree the proper- 
ties (i) of D, and will furthermore have good locality properties. To see this, 
note that if D is an exact solution of the GWR then the overlap formula just 
gives D back again: Dov = D [TB]. It is known that sensible ultra- local lat- 
tice Dirac operators cannot exactly satisfy the GWR jOji but approximate 
solutions are possible, and for these we have Dov ~ -D, indicating that Dov 
approximately inherits the properties of D and is also likely to have good 
locality properties since D is ultra-local. 

A specific construction of a short range lattice Dirac D with the prop- 
erties (i) and (ii) above arises from the perfect action formalism. This for- 
malism produces, in principle, a lattice fermion action free of artifacts in the 
scaling behavior via the iteration of renormalization group transformations. 
Moreover, the corresponding lattice Dirac operator satisfies the GWR as well 
jS]. The construction of the perfect action can only be carried out explicitly 
for free and for perturbatively interacting fermions |2I]-^ Still, at the non- 
perturbative level one can construct approximations of the classically perfect 
action (or fixed point action) of asymptotically free theories. The scaling 

^In that case, the axial anomaly takes the correct form even at finite lattice spacing 

m 
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artifacts of the fixed point action tend to be very small — as a study in the 
2d 0(3) model revealed j23] — and the fixed point Dirac operator solves the 
GWR too [7j. Hence approximations, or truncations, of the fixed point Dirac 
operator are natural candidates for operators with the properties (i) and (ii) 
above.® In particular, a hypercubic approximation of the fixed point Dirac 
operator has been considered, and it was found in numerical studies that the 
use of this operator as input in the overlap formula can significantly improve 
the scaling, locality and convergence properties of the overlap Dirac operator 
fHl CHI 1201 • These properties have been demonstrated in the Schwinger model 
[TH] . and also in QCD they have a potential to compensate the additional 
complication in Dhf compared to Dw: for the HF described in Ref. |TH] the 
locality of Dov is improved by a factor of 2 in the exponential decay com- 
pared to the standard overlap fermion. ^ Also the convergence rate increases 
significantly. However, to establish the overlap HF operator as a viable Dirac 
operator for lattice QCD one also has to check the conceptual basis, in par- 
ticular whether the correct axial anomaly and index are reproduced, and this 
is the issue that we address in the present paper. Experiments with simpler, 
non-standard operators inserted in the overlap formula have been performed 
in Refs. ^27|j. All those formulations are also covered as special cases by the 
considerations in this paper. 

The organization of the paper is as follows. In Section 2 we discuss the 
properties of D^f which are needed to compute its axial anomaly and index 
in Section 3. The conclusions and an outlook on further generalizations are 
given in the Section 4. 

^For direct QCD applications of a truncated fixed point Dirac operator, see Ref. j2()ll24j . 
However, the truncation distorts the chiral symmetry; indeed, even for truncated perfect 
actions the additive mass renormahzation can be considerable |25| . Chiral symmetry can 
be re-imposed though by inserting the truncated fixed point operator into the overlap 
formula, and this is another motivation to study the overlap operator with truncated fixed 
point operator as kernel. 

^This number refers to quenched QCD at /3 = 6, and the corresponding test for the 
standard overlap fermion was performed in Ref. |2(i) . 
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2 The structure of the hypercube Dirac op- 
erator 



We are going to use the following conventions for the 7 matrices: (7^)^ = 

7m, = '^S^iy, l2n+l = ■ ■ ■ l2n , SO that (72n+l)'^ = l2n+l and 

tr(72n+i7i ■ ■ ■ 72n) = (— 2i)". Throughout this paper summation over re- 
peated indices is implied. 

What we consider here is the minimally gauged HF-Dirac operator, which 
we are going to describe now. For techniques to simulate such HFs in QCD, 
see Ref. jigg . 

Assume the two lattice sites x and y to belong to the same lattice hyper- 
cube, i.e. — y^l < a for /i = 1 . . .2n. Then we denote by P{x,y) the set 
of lattice paths of minimal length connecting x and y. All these paths are 
inside the same hypercube again, their length is n^y G {0,1,..., 2n}, and the 
number of such paths is n^.^!. Let a be one of these paths. In the presence 
of a compact lattice gauge field we denote by U{a) the product of link 

variables along the path a. We also define the sign function e{t) = signt for 
t ^ 0, and e(0) = 0. Then the minimally gauged HF operator can be written 
as 

Dhf = ^(7mPm + A) , (2.1) 
= A„.,-^ E U{a). (2.2) 

ll-xy aeP(x,y) 

The parameters ki, . . . and Aq, Ai, . . . , A2n are coupling constants. Here 
we stay with the simple 7-structure of the Wilson fermion. ^ The vector 
term alone characterizes a generalized naive fermion, which would by itself 
generate species doublers. The scalar term A can be considered a generalized 
Wilson term which removes these doublers (resp. attaches a mass of the cutoff 
scale to them) for suitable couplings, see below. Hn^y and Xn^^ couple one 
site X to 2"^2' sites y. ^ 

®HFs with a more general Clifford algebra have also been used in QCD simulation, see 

Refs. EHEllEni- 

^The usual Wilson-Dirac operator with Wilson parameter r and bare mass too is re- 
covered by setting ki = 1/2 , Ao = 2nr+amQ , Ai = — r/2 and Kj = Xj = for j > 2. 
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Dhp has the correct formal continuum hmit with vanishing bare mass 
precisely when 

2n 2n 

$:2^-[7_-/]k, = 1, E2^T;]^. = 0. (2.3) 

At a finite 13 one deviates from these constraints and amplifies each coupling 
(except for Aq) in order to compensate its suppression by the hnk variable 
[T^ . However, in the current context we do impose the above constraints 
because they have to be restored in the classical continuum limit. They can 
be used to eliminate ki and Aq ; then Dhf contains 4n — 1 free parameters 

1^2) ■ ■ ■ 1 l^2n Ai, . . . , A2n- 

For our purposes it is useful to express Dhf in a coordinate-free way as 
follows. We use the parallel transporters 

T+^(x, y) = Uf,{x)5^^y_afL , T_f,{x, y) = Ul{x - il)6^_afi,y , (2.4) 
(where fi is the unit vector in /i direction) to define the Hermitian operators 

S, = ^^iT+,-T^,) , C, = ^iT^, + T^,) . (2.5) 

Moreover we introduce the following notation for the symmetrized product 
of operators, 

[01 • • • 0p\sym = - E • • • ' (2-6) 

P' a 

where the sum runs over all permutations a of {1, . . . ,p}. Then eqs. ()2.2j] 
can be re-expressed as 

2n 

~ i Pfi = — E 2''^p E [S^i'^ — C^^ ■ ■ ■ C^^)]sym , 

p=2 v2< - <yp ; t^j^tJ. Vj 

2n 

A = E2n-Ap) E [^-C.,---C.Xym- (2.7) 

Note that pj^ = — , A''" = A, hence Dhf satisfies 72n+i-Hermiticity, 

D'^hf = l2n+i Dhf 72n+i • (2.8) 

We now consider the zero-modes of the free field [U = 1) "naive" HF- 
Dirac operator ^'J^p^ and their "masses" provided by the scalar term ^A 



in eq. (|2.ip . The free field momentum representation of S*^ , (i.e. their 
eigenvalues for the plane wave eigenfunction e*'^^/") is obvious from eq. (j2.5p . 

Sf,{k) = sm{k^) = , C^{k) = cos{k^) = . (2.9) 

Hence the free field momentum representations of and A are 

2n 

p,ik) = ^^41-^2% (l-c,,---c,j] (2.10) 

p=2 y2<---<i'p; t^ji^iJ. Vj 

2n 

X{k) = En-Ap) E (l-c.,---c.J (2.11) 

p=l !/l<---<t'p 

The former vanishes when = 0, so the naive HF-Dirac operator has the 
usual zero-mode at = and the familiar 2^" — 1 "doubler" zero-modes 
for k at the corners of the Brillouin zone, just as in the case of the usual 
naive Dirac operator. However, in addition to these, there can be other zero- 
modes corresponding to vanishing of the factor in the square brackets on 
the right-hand side of eq. ()2.1U|) . It can vanish for some momenta k with 
components different from and vr, when K2, . . . , are in certain regions of 
the parameter space. These zero-modes correspond to new "exotic" spurious 
fermion species: if such a zero-mode occurs at k = k^^^ we set k = k^^^ + k' 
and find that the leading order term in the expansion of 7^p^(A;) around 

^(0) 

is ~ 7^ ^i^T^f^ ^'u- The corresponding propagator does not describe a 
usual Dirac fermion species. These exotic species are excluded though if the 
parameters K2, • • • , i^2n satisfy 

2n 

X/.N,...,K2n](/i;) = E2% E Cu2- ■ -Cvp) < 1 VA; . 

p=2 V2<---<^p; i^jT^M Vj 

(2.12) 

(Note that if this is satisfied for a particular index p then it is satisfied for 
all /i = 1, . . . , 2n. Also, since X/x = at /c = 0, the condition < 1 is the 
same as X/^ 7^ 1-) 

At this point it is natural to ask what are the values of ki, . . . , K2n that are 
of interest in practice, and do they satisfy ()2.12p ? The values of Ai, . . . , A2n 
used in practice are also relevant here since they determine \{k) and hence 

^"^Note that k represents a re-scaled momentum, which is 27r periodic at any lattice 
spacing a. Nevertheless we denote k simply as "the momentum" . 
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the masses of the doubler fermion species. As discussed in the introduction, 
one of the main aims in choosing the couphng parameters is to make Dh f as 
close as possible to satisfying the GW relation. The procedure of truncating 
perfect fermions (described in Section 1) led to the following values for the 
couplings |25| ITHj in dimensions 2n = 2 and 2n = 4:^^ 
2n = 2: 

Ki = 0.309 , K2 = 0.095 
Ao = 1.490 , Ai = -0.245 , Aa = -0.128 (2.13) 

2n = 4: 

Ki = 0.137 , = 0.032 , /t3 = 0.011 , /t4 = 0.005 
Ao = 1.853, Ai = -0.061, A2 = -0.030 , A3 = -0.016 , A4 = -0.008 

(2.14) 

In two dimensions the left-hand side of eq. (j2.12p is 

;^(2"=2)=4K2(1-C.)e^, . (2.15) 

The maximum of this, attained at Cu = —1, is 8k,2 = 0.76 for the coupling 
values in eq. (j2.13|) . hence eq. (j2.12j) is satisfied. In the dimension four case 
the left-hand side of eq. ()2.12|) can be re-written as 

r2n=4) _ i2K2 + 24K3 + 16ft4 + 16^(a3- (a + c^J(a + c^3)(a + c^J 

-16(^4 -)c^^c^.^c^^ , (2.16) 

where {^^,1^2,^3,^4} = {1,2,3,4} and a = ^2/(21^3). From this we see that 
when K2 > 2% (i.e. a > 1) and K4 > k\/ K2 the maximum of x^^"^^"* is 
attained at c^^ =c^^=c^^ = —1. For the coupling values ()2.14|) this maximum 
is 0.93, hence ()2.12|1 is again satisfied. 

Let {k^^^} denote the momenta of the zero-modes for the free field naive 
HF-Dirac operator. The mass of such a mode, provided by the scalar term 
in Dhf, is M^^^a where M^-^^ = X{k^^^). To avoid species doubling in 
the full HF-Dirac operator we impose the requirement on Ai, . . . , A2n that 

11 We give the values to 3 decimal places; they are given to higher precision in |25l I16| . 
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A(A;(^)) > for k^^^ ^ 0. From eq. (ETT|) we see that a sufficient condition for 
this is Ai < and Ap < Vp = 2, . . . , 2n, which holds for the couphng values 

For 
= 



in eqs. fl2.13|) - ()2.14|) (and of course also for the Wilson-Dirac operator) 
the usual zero- and doubler modes, characterized by s^ = V/i, i.e. 



or 71 for each fi , let N-,^ denote the number of fc^'s equal to vr. Then, from 
(j2.1ip . the mass M/a of the mode is seen to depend only on A^^ as 



2n 



E2--^^(-Ap)E 

p=l q=l 



L P- 









-1 . 







(2.17) 



with " =0 for q > Nj^. From this the masses in dimensions 2 and 4, 
with couplings given by eqs. (j2.13p - (j2.14p . can be determined (see also Fig- 
ures 1 and 2 in Ref. ^Hl)- We list them in the Table below in lattice units. 
For comparison we also list the masses of the corresponding modes of the 
Wilson-Dirac operator with Wilson parameter r. 



= 2: 





1 


2 






Mhf 





2.004 


1.960 






Mw 





2r 


4r 






= 4: 

iV. 





1 


2 


3 


4 


Mhf 





1.988 


1.960 


1.964 


2.000 


Mw 





2r 


4r 


6r 


8r 



The HF doubler masses are all close to 2 in lattice units, reflecting the 
fact that the free field Dhf with the coupling values from eq. ()2.13|) resp. 
fl2.14|) are good approximate solutions to the GW relation (since for exact 
GW solutions the eigenvalues lie on the circle in the complex plane centered 
at (1/a, 0) with radius 1/a). We also remark that, unlike in the Wilson case, 
Mhf does not always increase with increasing N^^. 

3 The continuum limit of the axial anomaly 
for the overlap-HF Dirac operator 

For a given value of the parameter m, the momenta of the zero-modes of the 
free field D^v are the k^^^ with M^^^ < m (both defined in Section 2). Hence 
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the parameter region in which Dg^ has a physical zero-mode and no doublers 
is < m < min{M(^) ^ 0}. 

Our aim now is to evaluate the classical continuum limit of the topologi- 
cal charge density q{x) given by eq. ()1.4|) . or equivalently, the axial anomaly 
A{x) = 2iq{x), of the overlap-HF Dirac operator. Specifically, we consider 
the situation where Dov is coupled to the lattice transcript of a smooth con- 
tinuum gauge field A = A^{x)dx'^, i.e. the link variables are 



U^{x) = Texp (^a / + (1 — t)afl) dtj , 

<^ 

where T implies t-ordering. We will derive the following result: 



(3.1) 



//m ^ M^-') Vj , then 



lim q{x) = I{k2, . . . , K2„; Ai 

a— >0 



where 



1 



Qcont y-^) 



(27ri)"n! 



^/^l.../^2r 



tr 



>^2n;rn) qcont{x) 



^X) ■ ■ ■ -^/i2n-lM2n (•^) 



(3.2) 



(3.3) 



is the continuum topological charge density and I{k2, ■ ■ ■ , ^2n', fn) is the degree 
of a certain map B : T^" S"^" given in eq. i3. 2^) below. In particular, 
I{k2, ■ ■ ■ , ^2n', ^) = 1 holds for m in the physical (doubler-free) region < 
m < min{M^^^ ^ 0}. Thus, form in this region, q{x) and the axial anomaly 
reduce to the correct continuum expressions in the classical continuum limit. 
Furthermore, when the parameters k,2, ■ ■ ■ , i^2n satisfy the constraint 
then 



/(k2, . . . , K2n; Ai, . . . , A2„; m) = 

{N^:M{N^)<m} 



2n 



(3.4) 



We first derive the result in the infinite volume, i.e. on a hypercubic 
lattice on IR^", and thereafter we discuss the finite volume 2n-torus case. 
The expression fll.4|) can be re-written as 



q{x) = -^tr(-^)(x,a;) 



iJ2 

m 



where 



Hm = l2n+l{aDHF - m) = -f2n+lilf,P^L + A - m) 



(3.5) 



(3.6) 



12 



is the Hermitian HF-Dirac operator (normalized by 1/a). We proceed as in 
the Wilson case treatment of Refs. [inillll by expanding (if^)~^/^ as a power 
series. First, H"^ is decomposed as 



Hi = L-V 

L = -p^ + {\-mf , 1/ = 7^[p^,A] + |7^7^[p^,p^] . (3.7) 

As in the Wilson case we observe \\V\\ ~ 0{a?), which is a consequence of the 
property || [T±^ ,T±y\ \\ ~ 0{a?) and eqs. (j2.7|) . Furthermore, a lower bound 
< 6 < L exists when the lattice is sufficiently fine. This implies that 
\\L-^V\\ ~ 0(a2) and consequently (Hl)-^/^ = (L[l - L-'^V])-^/^ can be 
expanded as a power series in L~^V when the lattice spacing a is sufficiently 
small. This was done in the Wilson case IT^ using an integral represen- 
tation of the inverse square root. (Note that the integral representation is 
needed since L and V do in general not commute.) The argument relies only 
on general properties of L and V which continue to hold in the present HF 
case; the treatment in Refs. ^3 generalizes straightforwardly to the HF 
and to arbitrary even dimension. Substituting the resulting expansion into 
eq. (|3.5|) and using the lattice (5-function 

^ ^"'''^^ ' My) = e^'^/" , (3.8) 

to express q{x) as an integral over momentum space, one obtains 

1 1 „7r ^2nr, i-yf fj (l,\p-ikx/a-[/n ikx/a\ 

' 'i^.r (3.9) 

where Hm{k) and L{k) are the free field momentum representations of Hm 
and L, and 

I dJ" f N-i/2 (2n)! 
c{n) = ——{l~t] = J , . 3.10 



^^This was established in the Wilson case for restricted values of m in Ref. [^ISOli and 
later for general m in j81| . The result will be generalized to the present HF case, and more 
general cases, in p?^ . 

^•^To derive eq. H;-{.9|l we have used the fact that the trace of the product of 72^+1 with the 
product of less than 2n 7 matrices vanishes. The factor in eq. (|3.9|) originates from 

the a^" in the operator representation 0iIj{x) = a^"-^y0{x,y)'4){y) (0 = Hm/ yjH'^ ). 
Hence the first term is of 0(1). 
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To evaluate the limit a — of eq. (|3.9|) we start from the following general 
observations, 

^-ikx/a [y^^ ^ y^^i giWa _ -a2F^^(x) e*^^'^'*^'^'^) + 0(a3) , (3.11) 
which imply 



g-,fcx/a[^^^^^]g.Wa = _a^F,,(^^)s^s, + 0{a') , (3.12) 

with the terms defined in eq. ()2.9|) . In the following we denote the free field 
momentum representation of a general lattice operator X by X{k). Then 
the relations (j3.12|) can be expressed collectively as 

[X , F] e*'^"/" = -a^F^(3ix)daX{k)dpY{k) + 0{a^) (3.13) 

for X = and Y = Si, , C^. In fact this relation continues to hold when 

X and Y are general polynomials of the 5*^ and the C^. Since and A are 
such polynomials we can apply eq. (1^11^ to e"^^^/'^ 1/" e*^^'/'' in the expression 
fl3.9p to obtain 

= i a C/ii.../i2n ^^(-^aia2('^) ■ ■ ■ -^a2n-ia!2n('^)) ^ 

(A(A;) - m) 9„,p^,(A;) ■ ■ ■ d^.,^p^^Sk) 

-2n9„,p^,(A;)---9«,„_,p^,„_,(/.)9„,„A(/.)p^,„(/.)] +0(a2"+i) . (3.14) 

We now note the two identities, which will be crucial for our further consid- 
erations: 

eMi.../^2„ tr(F^,^2(x) ■ ■ •F/.2„„i/.2„(x)) e„i,...«2„ (9„iPi(/i;) ■ ■■da2„P2n{k) , (3.15) 

C/ll.../t2n tr^^-F^jaj (x) ■ ■ ■ -fa2ri-ia2n (''')) ^ 

daiPf,^{k) ■ ■ ■9a2„_iPM2„-i(^) 5a2„A(A;)p^2„(A;) = 

2n 

p=l 

daoKk) da^pi{k) ■ ■ ■ 9«j,_iPp_i(A;) 9q^+iPp+i(/c) ■ ■ ■ 5„2„p2„(/c) (3.16) 

14 



These combinatorial identities rely only on the facts that the p^{k) and X{k) 
all commute, on F^y{x) = — F^^(x), and on the cyclic property of the trace. 
Replacing the left-hand sides of these identities by the right-hand sides in eq. 
()3.14|) . and introducing the real- valued functions 

eQ{k) = -[X{k) - m] , e^{k) = -tp^{k) (/i = l,...,2n), (3.17) 

we arrive at 



tr 



2n 



J2{-lY OpdOoA-.-A dOp^i A dOp+i A ... A d^2n + 0(0^"+^) , (3.18) 

where d6j = daOjdka is the exterior derivative of 6j. Substituting this into 
eq. ()3.9|) and re-writing L from eq. ()3.7|) as 

2n 

m = -p\k) + [A(A;) - m]2 = 5: O^ikf = \e{k)\^ , (3.19) 

p=0 

we obtain hm^^o li.^) = I{.i^2, ■ ■ ■ , A2n; rn) qcont{x), as we claimed in eq. ()3.2|) . 
with 

1 n' 
/(k2, . . . , A2„; m) = - c{n) — X 

1 ^" 

/ ^^E(-l)''^P^^oA... Ac/Vi Af/^p+i A... Arf^2n . (3.20) 

A little calculation shows that the integrand here can be re-written as 

£(-1)4 4^) A ... A A A ... A . (3.21) 



p=0 



This is precisely the pull-back to T^" =] — vr , vr]^" of the volume form on the 
unit 2n-sphere S"^" C |R^"+^ via the map 
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Furthermore, the coefficient of the integral in eq. (|3.2U|) turns out to be (recall 
definition ^TW^ ) 

where ^0/(5*^") is the volume of the unit 2n-sphere. Hence expression (I3.2(jj) 
calculates the degree of the map (j3.22|) . 

I{K2,...,X2n,m) = deg{e) . (3.24) 

This is a generalization of the topological evaluation of the anomaly coeffi- 
cient given in Ref. 



If ||| G S'^"' C |R^"+^ is a regular point for the map B then it is a standard 
topological fact that deg{Q) = J2iSi , where s/ = ±1 is the sign (relative to 
dki A ... A dk2n) of the integrand of eq. (j3.2Up evaluated at a pre-image 
of -jli , and the sum is over all the pre-images (labelled by /). We choose 

1^ = (1, 0, . . . , 0). Then the pre-images are precisely the subset of the 

k^^^ introduced in Section 2 which satisfy 9o = — [\{k) — m] > 0, i.e. for which 
M^^^ = \{k^^^) < m. Moreover, the integrand in eq. ()3.20|) reduces at these 
momenta to dOi A ... A d92n- To determine the sign of this at a given k^^\ 
recall from eq. ()2.1()j) that 9fj,{k) = —ip^{k) = S^{k)[l — with x^l^) 

defined in eq. fl2.12|l . It follows that 

2n 2n 

de,A...Ad92n = {l[C,ik))[l[il-x.ik)))dk,A...Adk2n 

+ terms with at least one S^{k) factor. (3.25) 

If < m < min{M'^''^ ^ 0} then there is precisely one fc*^-'-' for which M*^-'^ < 
m, namely fc*^-^^ = 0. In this case, since — ^ ^ fj. , eq. (j3.25p yields 

d9i A ... A d92n = dki A ... A dk2n , i-e. the sign is +1, so deg{Q) = 1 in 



k=0 



k=0 



this case, as we claimed. Let us now consider the case where the parameters 
K-2, ■ ■ ■ , i^2n satisfy the constraint ()2.12|) . i.e. Xfii^) < 1 W fj, ,k. Then the k^^^ 
are precisely the k at which Sf^{k) =0 V/i, so the terms with S^{k) factors 
in eq. (j3.25j) vanish and the sign of the remainder is given by n^=i Cfj_{k^^^). 

This sign is (—1)^^^' where Nj^^^ is the number of components of fc^-^-* which 
are equal to vr. Recalling from Section 2 that M^^^ depends only on iV^-^^ 
in this case, and noting that the number of k^^^ with A^^-'^ = A^^^ is 



2n 
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it follows that deg{Q) = J2{n^ ■.M{N^)<m} n {~^) ■ This completes the 
derivation of the claimed result ()3.2|) - ()3.4j) in the infinite lattice setting. 

The rigorous derivation of the expansion ()3.9p . carried out along the same 
lines as in the overlap- Wilson case J3 E], requires an assumption on the 
continuum field, namely that Af^{x) and its first few partial derivatives are 
bounded on IR^". Such bounds are guaranteed to exist if A has a compact 
support on IR^". However, having established the result in eqs. ()3.2|) - ()3.4|) 
for gauge fields with compact support, it can then be extended to general 
smooth gauge fields using locality type arguments, in the same way as in the 
overlap- Wilson case (see eq. (3.45) in Ref. [T3j and the associated discussion). 
This relies on the existence of a non-zero lower bound on . 

In the finite volume 2n-torus setting the momentum integrals in eqs. 
fl3.8|) . ()3.9|) become sums, so the derivation given above does not carry over 
directly to that setting. However, using locahty-based arguments one can 
show that the finite volume q{x) coincides with the infinite volume q{x) 
up to exponentially suppressed finite size effects, thereby establishing that 
the continuum limit results ()3.2j) - ()3.4|) continue to hold in the finite volume 
setting. This was done in the overlap- Wilson case in Ref. jT3] ; the arguments 
there relied only on general properties and carry over to the present HF case 
(given the aforementioned lower bound on H^); the details of all this will be 
given in the general setting in Ref. |32j . 

In the light of the index formula (jl.7|) it then follows that when the 
overlap-HF operator is coupled to the lattice transcript of a smooth contin- 
uum gauge field A on T^" with topological charge Q, then index{Dov) reduces 
to I{k2i ■ ■ ■ 1 ^2n'i fT^) Q in the limit a — > 0. Thus the HF fermionic topological 
charge reduces to the continuum topological charge in the classical contin- 
uum limit when the parameter m is in the doubler-free region, just as it does 
in the Wilson case. 



4 Summary 

We repeat that the previous literature contains the following considerations 
about the axial anomaly of overlap fermions: 

• Perturbative considerations show that the correct continuum limit is 
obtained in the sector of topological charge zero for any overlap oper- 
ator, see in particular Refs. [TUj . 
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• There was also a rigorous, non-perturbative proof that covers all topo- 
logical sectors, but it was specific to the case of the simplest stan- 
dard overlap fermion, which uses the Wilson-Dirac operator as an input 

The standard overlap operator is wide-spread in recent simulations. How- 
ever, there are attempts by various groups to use also non-standard overlap 
operators [IHl dl 1201 12Z] • The operators used in those works are all included 
in the class of HF overlap operators. For the latter we have given in this paper 
a non-perturbative evaluation of the continuum limit of the axial anomaly 
and index which is valid in all topological sectors. 

We have formulated the HF-Dirac operator in 2?2-dimensional Euclidean 
space in the form ()2.7|) . which is well-suited to analytic investigations. We 
used it first to study the dependence of the doubler structure of Dup on 
its coupling parameters. Then we evaluated the classical continuum limits 
of the axial anomaly and index of the overlap-HF Dirac operator, show- 
ing that the correct continuum expressions are recovered when parameters 
are in the physical (doubler-free) region. A noteworthy feature of our con- 
tinuum limit evaluation is that it relies only on general properties of the 
HF-Dirac operator and not its explicit form. This is in contrast to the pre- 
vious evaluations in the Wilson case (a special case of the more general 
HF structure) which all use the explicit form of the operator. The main 
new technical observations which our approach is based on are the general 
relation ()3.13p and the identities ()3.15|) - ()3.16|) . These ingredients allow the 
continuum form e^j...^2^ tiF^^^^{x) ■ ■ ■ F^,^^_^^^^{x) of the axial anomaly to be 
extracted, and its coefficient to be topologically evaluated as the degree of a 
map : T^" 5*^" using only general properties of the HF-Dirac operator. 
These properties are not specific to the HF case, and the approach can be 
extended to completely general overlap Dirac operators obtained by substi- 
tuting a general ultra-local lattice Dirac operator (involving the full Clifford 
algebra of 7 matrices) into the overlap formula (jl.8|) , as it was done in Ref . 
pUj . The full extension, which involves considerable additional work, will be 
carried out in a forthcoming paper |32j, and it will be shown there that the 
axial anomaly and index of the general overlap Dirac operator Dov continue 
to have the correct classical continuum limits in the physical (doubler-free) 
parameter region specified by < m < minjM^-'^ 7^ 0}. It was worthwhile 
to consider the overlap-HF case on its own, firstly because of the current in- 
terest in using this operator in numerical simulations, and secondly because 
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it illustrates the main ideas and techniques of the general case but without 
the extensive formalism and additional complications of the latter. 
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